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0. ABSTRACT 


est linear unbiased estimation of location and scale parameters of the 
Weibull distributions using ordered observations of a random sample is con- 
sidered. It is assumed that the shape parameter of the Weibull distribution is 
known. For sample sizes up to and including five, all possible censoring is 
considered. For sample sizes greater than five, one-sided censoring (that is, 
large values of the sample are missing) is considered. The coefficients are 

fabsar fern: 
presented in . For each sample size and the value of the shape para- 
ha Gee 
meter, the first row of coefficients in fabite I correspond to the best linear un- 
a 

biased of the location parameter and the second row of coefficients corres- 
pond to the best estimation of the scale parameter. The accuracy of these co- 
efficients is to four or more decimal places for sample sizes less than or 
equal to 5, to three or more decimal places for sample sizes 6 to 9 and to two 


_— 


or more decimal places for sample sizes 10 to 12. ee S@Z-nas also com- 


puted these coefficients for sample sizes up to and including 15, using ap- 
proximations for the variances and covariances of the Weibull order statistics. 
However, for sample sizes greater than 2, our values for the coefficients do 

< not agree with his values. This is not surprising since the approximations 


for the variances and the covariances used by Weibull [3€}-are not very ac- 


curate, 
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The expected values and the variances of the Weibull order statistics 
are computed accurate to four or five decimal places. These are presented in 
Table 2. In Table 2, for each combination of N and i the first row gives 
the expected values and the second row gives the variances. The product moments 
and the covariances are presented in Table 3. The accuracy is four to five 
decimal places except for sample sizes 11 and 12 where the accuracy is three 
decimal places. Again, for each combination of N, i and j, the first row 
gives the product moments and the second row gives the covariances. Each 
column of values in Tables 2 and 3 corresponds to a certain value of the shape 
parameter. Throughout, the same values of practical interest are considered 


for the shape parameter of the Weibull distribution. 
1. INTRODUCTION 


Many contributions have been made to the problem of linear unbiased es- 
timation of the parameters of a distribution using the ordered observations of 
the sample. Lloyd [8] tackled this problem in general and gave explicit form- 
lae for estimating the location and scale parameters of a distribution using 
order statistics in a random sample. Gupta [3] considered the problem of best 
linear estimation of the mean and standard deviation of a normal population using 
a censored sample. Sarhan [10, 11, 12] considered the problem of best linear 
estimation of the mean and standard deviation of a population by order statistics. 
Sarhan and Greenberg [13, 14] tackled the problem of estimating the location and 
scale parameters by order statistics in singly and doubly censored samples from 
normal and exponential populations. Tate [16] considered the best linear unbiased 
estimation of functions of location and scale parameters of a distribution. 
Plackett [9] gave some approximate expressions for the variance-covariance 


matrix of the censored sample and thus gave some 'almost' best linear estimated 


of the parameters of a population. Blom [1] gave some approximations for the 
expected values, variances and covariances of order statistics in samples 
drawn from a continuous population such that the inverse of its cumulative 
distribution function could be expanded in Taylor series. Gupta [4] gave the 
moments of the order statistics in a random sample drawn from a gamma popula- 
tion with integral values for the shape parameter. He also considered the 
best linear estimation of the location and scale parameters of the gamma popu- 
lation. Weibull [17] gave a distribution which has been widely used in life 
testing problems. Lieblein [7] gave exact closed form expressions for the 
yeh moments and the product moments of the Weibull order statistics. Weibull 
[18] by applying the approximate formulae of Blom [1] for the variances and 
covariances of Weibull order statistics, computed the coefficients for the 
best linear estimation of the location and scale parameter of the Weibull dis- 
tribution using the complete sample. However, the coefficients are in great 
error, Leone et al [6] considered method of moments and method of maximum 
likelihood for estimating the parameters of the Weibull. Dubey [2] considered 
approximations to the maximum likelihood estimators of the parameters of the 
Weibull distribution. White [19] recently showed that by considering the 
logarithms of the Weibull order statistics one can obtain best linear unbiased 
estimates of the scale and shape parameters of the Weibull distribution as- 
suming that the location parameter is zero or known. Numerical values of the 
moments of Weibull order statistics have not been tabulated so far. Hence, it 
is of interest to numerically compute the lower moments of Weibull order sta- 
tistics for small sample sizes at least and obtain accurate and best linear 
coefficients for the estimation of the location and scale parameters of the 


Weibull distribution. 
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3. NOTATION AND KNOWN RESULTS 


In this section we will give the least square estimates of the loca- 
tion and scale parameters of an arbitrary continuous distribution as given 
by Lloyd [8]. However, we will consider general formulae that take care of 
any kind of censoring. We also give the exact formulae of Leiblein [7] for 
the moments of the Weibull order statistics. 

Let X be a random variable with the density function f((x-6)/¢) so 
that the transformed variable Y = (X-0)/a has density f(y). That is, Y 


has Oandl1 for its location and scale parameters respectively. Let 


Xx <x < vee CK, 


ry ,N 5 “r,+1,N 5 l¢r 


aN? 1S <0, (3.1) 


be the available portion of the ordered sample of size N drawn from the 


population with density f((x-6)/s). That is, the first r,-l and the last 


1 
N-r, observations are either missing or ignored. We wish to obtain the 
best linear unbiased estimates of © and «o on the basis of this censored 


sample. Let 
Ys on = (Xs iy - 0)/s, A= ry, ry +1, sees Py, (3.2) 


and 


Man Bs WS 5,4 n = COVEN ye Eye ESAS dS FQ (3.3) 
We assume that Us oN and 4 4,N are known. Then 


= _ 2 
E(K, yy) =@+ 6 sy and Cov(X, vy» Xs ow) =6¢ 34 5,N 7 


' Since E(X, n) is linear in @ and og, by Gauss-Markov theorem, @ ando are 
? 


estimable unbiasedly by linear combination of X, Let 
$ 


N° 


1 


= ' = ! 
xX mA GoW? eeey Xow ’ I oe (1, @eey 1) 


= § s 


t -t 
Ww = (up uN? eeey Mr 5 and Q = (54 5 5) bd 


Then, in matrix form, we have 


E(X) = Ile +o yw = ((I,y)) ) 


and 2 (3.4) 
Var(X) = 6“w, 


where w = qr 


is an (ro-r, +1) x (ro-r,+ 1) positive definite symmetric 
_matrix. Denoting ((I,y)) by P and (@,6) by &', we have by Gauss-Markov 


theorem (See Lloyd [8]) the estimate of «; 

& = (p'a p)7) p'a x = ox (3.5), 
where bisa 2x (ro-r, +1) matrix which will be called the coefficient matrix. 
Also, the variance-covariance matrix of the estimates is given by 


(P'a p)7+ 6? , (3.6) 


where : 
POT. Dog 


May u'@u 


P'Q P= 


the elements of the matrix being scalars, The inverse is given by 


t 
(P'Q p)~+ = avi wpQuy ae Th 
-I'Q y r'Qi 
where A is the determinant of the matrix P'QP, Using this result in (3.5) 
and (3.6) we obtain, 


@=-y'DX, $ =1'DK (3.7) 
where D is the skew symmetric matrix defined by 


D=Q(I y'= yI')2/a (3.8) 
Also 
Var(§) = i Gia “JA; var($) = I'ets?/A 


and (3.9) 
Cov(8,$) = - I'ays?/a 


All the preceding formulae are explicitly given by Lloyd [8]. 
Next we consider the general Weibull distribution and assume that 
its shape parameter denoted by m is known. Consider the following distri- 


bution function for X. 
F(x) = 1 - exp(- {(x-@) /s}”), m>l, x>e. 


Weibull [18] remarks that for m< 1, the Weibull distribution has no practical 
value since the derivative of its density function is infinitely large at 
x = @. Thus, throughout,we consider only values greater than unity for nm. 


If Y = (X-9)/s, then the distribution function of Y denoted by G(y) is 


G(y) = 1 - exp(- y”) y>0O. 
Let Yh < YON See < Yy yu be the order statistics in a sample from G(y). 
Lieblein [7] has given explicit formulae for the x” moment of the Ys yy and 
’ 
the product moments of Yy ys and Tsou’ Thus 


i-1 
BYE y) = ce) r (atm) (02) (ietehay Hea) 
bro 


l<ic<N, k=1, 2, .. and m> 0, (3,10) 


and 
i-1 j-i-1 


BY Wy 150) =e i , (yb cy) (ty (j-it, -s,N-jtstl), 
Ls 8=0 L 


l<i<j<n, m>0 (3.11) 


where K = Nip(2+2m)/(4-1)4(j-4-1) 8(N-j) 4, 
=i 

y(tsu) = (ta)"*™ BB (lem, lim”), 
1 1 
=i. “i Pon n 
p=t/(ttu) and 8, (1m » ltm’) = ij x (1-x)" dx. 
0 


(3.10) is slightly different from the expression due to Lieblein [7]. However, 
this can be obtained by expanding gt-hy) in powers of 1-G(y) and integrat- 
ing termwise. The covariance between Y, N and a y can be obtained by 

> ’ 


subtracting E(Y, Welz from E(Y However, when m=1, the 
’ 


30) iN %,n)° 
Weibull distribution becomes the negative e.ponential distribution in which case, 


it is known (see Sarhan [10]) that 


E(Y 


i 
aos ), (een)? 


B=l 
(3.12) 


4 
Var(¥, y) = ), (wef)? 
=] 


and 


Cov(Y, = var(Y, yy) ; i<j. 


ww? yy) 


Equation (3.12) could be used for checking the accuracy of the moments of the 


Weibull order statistics. It is also well known that 
N 
) e(Yy y) = NE(YE) . (3.13) 
’ 


Consequently N 
), By y) = 8 rQ/m) 
i1=1 


and, 
N 


), BCE y) = N r(2/m) 
i=l 


which could be used for checking the accuracy of the numerical values of the 


moments of Weibull order statistics when m/l. 


4. TABLES PREPARED AND EXPLANATIONS 


Table 1: The best linear coefficients for the ordered observations 
X, N from the Weibull distribution given by 
’ 


F(x) = 1 - exp(- {(x-@)/o}”) 
are presented in Table 1 for m= 1.5(0.5)3(1)8 and for sample sizes 
N = 2(1)12. For each combination of N and m_ the first row of coefficients 
corresponds to the estimation of @ and the second row corresponds to the 
estimation of o. An easy check of these coefficients is that the coeffi- 
cients for estimating 6 add up to 1 and the coefficients for estimating o 
add up to 0. The accuracy of the coefficients presented in Table 1 is as 


follows: For N = 2(1)5, all possible values for r, and r, are con- 


sidered and the accuracy is four or more decimal places. For N = 6(1)9, 
r= 1 and only one-sided fixed percentage truncation is considered. Here, 
the accuracy is 3 or more decimal places. For N = 10(1)12, r= 1, and 


the accuracy is 2 or more decimal places. 


Table 2: The expected values and the variances of the Weibull order 
statistics in samples of sizes up to and including 12 are presented in 

Table 2 for m = 1(0.5)3(1)8. For each combination of N,i and m, the 
first entry gives the expected value and the second entry gives the variance. 
The accuracy is to 5 decimal places for N = 2(1)9 and to 4 decimal places 
for N = 10(1)12. r(1+#km"+) was computed following a procedure suggested 


by Haynam and Hansen [5]. 


Table 3. The values of E(Y, oy Ys ow and Cov(Ys ys Y5 Nn) for sample 
sizes up to and including 12 and for m = 1(0.5)3(1)8 are presented in 
Table 3. Again, for each combination of N, i, j and m, the first entry 


gives E(Y, , Y, ,) and the second entry gives Cov(Y The 
i,N j,N i 


ww? 45, 


accuracy is to 5 decimal places for N = 2(1)9, to 4 decimal places for N = 10 
and to 3 decimal places for N = 11 and 12, The following representation was 


used in order to compute B (wow) occurring in (3.11): 


- + 
Ba(wiw) = So (-2)® (Ht) e8/cwte) (3.14) 
p s 
s=0 
the summation is finite if w is a positive integer. When w is not an integer 
terms which are of magnitude less or equal to 1076 were neglected. 
Throughout, the following rounding off procedure was adopted. If 


-(nt1) was added to the 


accuracy is desired to n decimal places, 5.10 
computed value and n decimal places of the resultant value were retained 


for final presentation in Tables 1, 2 and 3. 


5. COMPARISON WITH WEIBULL'S TABLES 


Weibull [18] computed the expected values of the Weibull order 
statistics accurate to four decimal places using Leiblein's [7] formula, 
for m=1(1)10. However, there is a marked disagreement among his coefficients 
and our coefficients given in Table 1, especially when N is greater than 2. 
To get an idea of the disagreement, we will present the values for N = 5, 


m=5 inTasa-- 5.1. 


TABLE 5.1 
Weibull: for §: 1.79 0.67 0.31 -0.12 -1.65 
Our for 6: 1.56930 1.15448 0.33735 -0.62650 -1.43460 
Weibull: for 8: -1.69 -0.54 -0.15 0.33 2.05 
Our for 8: -1.69774 -0.96413 -0.08339 1.07556 1.50294 


ry nS Se 
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The discrepancy is not surprising at all, since Weibull [18] used crude 
approximations of Blom [1] for the variances and covariances of the Weibull 


order statistics. 


6. CONCLUSION 


The best linear coefficients and the moments of Weibull order 
statistics presented in Tables 1, 2 and 3 will be useful in life testing 
problems. All the computations were done on Burroughs 220 compiler using 
single precision arithmetic. The computation of these tables were completed 


in about nine hours of machine time. 
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